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Abstract 

We evaluate the long distance weak contribution to the neutron and 
proton electric dipole moments using an effective Lagrangian framework. 
We estimate the coefficients needed by a factorization hypothesis and ad- 
ditional assumptions on 75 terms in the baryon lagrangian. We obtain 
|dg I ~ 5 X 10"^'^ e-cm and \d^\ ~ 4 x 10~^^ e-cm. The former estimate is 
similar to the quark model estimates done previously. 
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1 Introduction 



The contributions to the electric dipole moment (EDM) of the neutron within 
the Standard model have been widely explored. In the Standard model the 
main contribution comes from the strong CP-violating ^-term [|[] and recently 
reconsidered in [^. The experimental and theoretical situation can be found 
reviewed in ||^. The present experimental limits are 

|ci^| < 11 ■ IQ-^^e ■ cm [4] and < = (-3.7 ± 6.3) ■ lO^^^e ■ cm § . (1) 

The neutron and proton EDM appear as an higher order effect of weak inter- 
actions. The one loop contribution with W exchange vanishes because of KM 
combinations and the two-W boson loops contribution is also vanishing as shown 
in and references therein. First non vanishing contributions are the so called 
transition quark electric dipole moments Q and the insertion of penguin 

diagrams within the baryon . Penguin diagrams can in fact produce the CP 
violating phases needed to generate the EDM term. The EDM is then gener- 
ated by a two step process: the strong penguin diagram insertion which causes 
the transition d — > s and weak radiative decay of the final strange baryon (e.g. 
I1^,A — > ri'-f). Already in [jlO| it was observed that penguin diagrams' contribu- 



tions dominate the EDM. The evaluation of the long distance part of penguin 



insertions has been done up to now relying on quark models, like the one in [ p!l[| . 
See for more references. 

In this letter we propose an alternative derivation of the neutron and proton 
EDM based on a factorization hypothesis which leads to the derivation of the 
EDM within the framework of chiral perturbation theory for baryons. 

We first describe our approach and perform the calculation. Here the assump- 
tions made at various stages will also be explained. Then we present numerical 
results for both the proton and neutron electric dipole moment. A compari- 
son with power counting in the heavy-baryon formalism for Chiral Perturbation 
Theory and a proof that our contribution is the leading one are presented next. 
Finally, we recapitulate our main conclusions. 

We do not attempt to ascribe an uncertainty to our results. However, contrary 
to the p-wave hyperon nonleptonic decays we have rather small cancellations 
between the different subamplitudes. We therefore expect higher orders to be 
of normal size. The main uncertainty is the assumption made in estimating the 
coefficients in the Lagrangian and the factorization ansatz. 



2 The calculation 

The gluonic penguin is the main source of CP violation in the weak \AS\ = 
1 hamiltonian. The effective interaction which mediates the d s transition 
in the EDM diagram involves the strong penguin four-quark operators of the 
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effective weak [AS"! = 1 Lagrangian £e// = —Gf/\/^ Ylt=^ Ci{jj)Qi{n) (we 
use tlie definitions of tlie Qi as in |12|)- Tlie operator Qq is defined as Qq = 
{sadi3)v-AJ2q{(li3(la)v+A, wliere V ±A stands for tlie combination 7^(1 ±75) and 
a, {3 denote colour indices. Using Fierz identities one can rewrite Qe as follows 

Qq = -^Y. So^^qan QPudh + 4^ Sa^l^^qau Qh^f^dfSL > (2) 

wliere qr^l = (1 ±75)/2 q. At this point we introduce a factorization hypothesis: 

(EilQel^,) = -8{dndL){B,\sLdn\B,) . (3) 

Within the factorization hypothesis all the other operators Qi, i = 3, 4, 5 do not 
contribute. The hypothesis is favoured by the substantial enhancement of the 
coefficient of the operator Qq by next-to-leading corrections and the enhancement 
of the Qq contribution to weak non leptonic decays. For later use we introduce 
qX±q with A-t projection matrices defined in terms of the Gell-Mann matrices 
Ag, A7 as A-t = (Ag ± 1X^)12. 

The baryon Lagrangian for strong interactions in presence of external scalar 
and pseudoscalar sources allows for the general form leading in the derivative 
expansion and in the light quark mass matrix: 



Cb = Tr[B{iD^r -mB)B) 

+hTr (Bx+b) + b^Tr (BBx+) + hTr {bb) Tr 

+h\Tr {BX-I5B) + blTr (Bj.Bx-) + blTr (Bj.b) Tr (x-) , (4) 

where Tr stands for the trace over flavour indices and the covariant derivative in 
the case of interest contains the electromagnetic field D^j^B = dfj^B + ieA^[Q, B]. 
The field x± = ^^X'C^i^X^C contains the external scalar and pseudoscalar sources 
with X = 2-Bo('S(a;) + ip{x)) and, in the absence of meson field = 1), we have 
X+ = 4:Bqs{x) and X- = 4:iBQp{x). Bq is related to the scalar quark condensate 
through the identity (0|gg|0),=,,,,, = -B^{p/2){1 + 0{M)), with / ~ ~ 132 
MeV and M. the light quark mass matrix. The baryon fields are incorporated in 
the 3x3 matrix 



B 



and transform non linearly under SU{'i)L x SU{?))b. The electromagnetic cou- 
pling of the baryon field in the covariant derivative generates the ordinary mag- 
netic moment for baryons. 
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Taking the functional derivative with respect to Sij,pij of the generating func- 
tional of the baryon Lagrangian including terms (^) it follows that 

{Bk\qiqj\Bi) = {Bk\ - ^BohB.Bj - ABob2{BB)ji - ABohTrBB5,j\Bi) (5) 
{Bu\q^-i^qj\Bi) = {Bk\AB^hlBa^B, + AB^hl{B^^B)^i + AB^hlTrB^^B5,,\Bi) . 

Both identities relate the two fermion matrix element @ to the corresponding 
baryon Lagrangian. In the case of the flavour combination BiBj with i,j = s,d 
also the ordinary baryon Lagrangian for \AS\ = 1 weak interactions has to be 
taken into account. Assuming octet enhancement it has two terms at leading 
order in the derivative expansion which transform as (8^, Ir): 

Clff^=' = a TrB{Xe, B} + b TrB[K, B], (6) 

where for the case of interest the meson field is absent, i.e. ^ = 1. Then the full 
baryon Lagrangian which induces s ^ d transitions can be written as follows 



C^<^ = 4Bo ( ^^^hTrBX_B + ^^^b2TrBBX. 
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+45o( ^—^hiTrB\+B + — —b2TrBBX^ 



-4B,^(blTrBX.^,B + hlTrB^.BX^^ 



+45oy [h\TrBX+-i^B + blTrB^,BX+) . (7) 

The parameters 7, 6 are defined as 2i?o7&i = a + b and 2i?o(5&2 = a — b, where a, b 
are the couplings in (^). The Lagrangian (0) induces transitions n S°,A and 
p S+ both parity conserving (p.c.) and parity violating (p. v.) i.e. with a 75 
insertion. Here we have made the assumption that the parity violating part of 
the other operators besides can be neglected. The EDM term of the neutron 
C-EDM = i^i'CTfj^^F^'^'y^'tp can be generated at tree level in the baryon theory by 
the vertices in (|^) and through the insertion of the anomalous magnetic moment 
operator which we write as 

^"^^ {^'DTrBa^,F^'''{Q, B] + AfipTrBa^^F^" [Q, B]^. (8) 

The baryon magnetic moments receive the ordinary contribution from the leading 
electromagnetic coupling in and the anomalous contribution from the next to 
leading O(p^) terms in 

At tree level in the full baryon theory the set of diagrams which contribute 
to the electric dipole moment of the neutron and proton are shown in Fig.|l[ 
They are given by the insertion of a parity violating (p. v.) vertex and a parity 
conserving (p.c.) vertex from and the insertion of the anomalous magnetic 
moment vertex in M)- 
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Figure 1: Tree level diagrams which contribute to the neutron electric dipole 
moment. The circle vertex is the parity conserving vertex. The box is the parity 
violating vertex. The photon insertion is the anomalous magnetic moment vertex. 
In the proton case the same diagrams contribute where the p S+ transition is 
allowed. 



The Lagrangians (|^ and (||) lead to the following tree level contribution for 
the neutron case: 



ie a — a* ( a a5 1 



2m AT 



1 



M, 



A 



^M^{Ml-Ml) 

5 Mso 

M„(M|o - Ml) 
1 



(Ma-M„)(Mso + M„ 
1 

(Mso -M„)(Ma + M„ 



(9) 



The last two terms are the mixed S° — A exchange contributions. The /Zj are 
the magnetic moments of the neutral baryons in units of nuclear magnetons (1 
nuclear magneton=e/2m7v = 1.052 ■ 10~^^ e-cm) and we assumed as valid their 
tree level expressions in terms of fiD'- 
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1 



(10) 
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while the parameters Ai short hand notation for 
As = -V2Bo6b2 A| = 2V2Bobl 
Aa = y|5o(5fe2 - 27&1) Al = -2^Bo{bl-2bl) 

(11) 

and are related to the weak non leptonic hyperon decay amplitudes mr" , 

A — > p7T~ as it is explained in the next section. Note that the anomalous magnetic 
moment term with coupling is the only contribution to the magnetic moments 
of neutral baryons. 

In the proton case the only possible transition induced by (0) is p ^ It 
involves only 62,^2 type of coefficients. In the charged case also the anomalous 
magnetic moment term proportional to A^p does contribute. The expression for 
the electric dipole moment of the proton induced by Lagrangians (0) and (H) is 

Here Afip = A/is+ = /Wp— 1 = l/Sfio+A^p are the anomalous magnetic moments 
of the proton and the S"*" in units of nuclear magnetons which are equal in the 
SU{3) limit. They receive contributions from A^p and fiu. The proton EDM 
involves only the anomalous magnetic moment terms. It is easy to verify that 
the tree level diagrams as in Fig.^ where the anomalous magnetic moment vertex 
is replaced by the ordinary electromagnetic coupling in the proton case sum to 
zero. 

The expressions and ([T^) are the tree level contributions to the electric 
dipole moment of the neutron and the proton respectively. In section ^ we ex- 
plicitly show the power counting for the tree level and quantum corrections in 
the Heavy Baryon Chpt and that one loop corrections are naturally suppressed 
also in virtue of the absence of large cancellations at tree level. 



3 Numerical results 

The numerical estimate of the tree level contributions to the neutron electric 
dipole moment in (^) and the proton one in requires the knowledge of the 
following set of weak and strong parameters: the combination a — a*, the pa- 
rameters v4a,s, ^as finally the magnetic moment coefficients fiD,Afj,F. 

The latter can be extracted at tree level from the measured values of the mag- 
netic moments of baryons. At tree level the SU{3) symmetric Coleman- Glashow 
relations amongst magnetic moments are valid, while they are experimentally vi- 
olated by about 0.25 nuclear magnetons in average [0. For the leading quantum 
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corrections to the magnetic moment of baryons in the Heavy Baryon expansion 
see e.g. 
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If we use the experimental values of /ip and yU„ to determine fxn and Afxp 
and disregard quantum corrections we obtain fiu = —\^n = 2.87 and Afip = 
jjLp + |/i„ — 1 = 0.8365 for the experimental values /Xp = 2.793, /i„ = —1.913 in 
units of nuclear magnetons. 

The other magnetic moments in the tree level approximation are: 

/iA = i/i„ (-0.613 ± 0.004 = -0.96) 
1 

^^20 = -^Ain (±1-61 ± 0.08 = 1.66) 

/i2+=/ip (2.458 ± 0.010 = 2.793). (13) 



In brackets the latest experimental values are indicated [|1J] and compared with 
the SU{3) symmetric value. The S° magnetic moment has not been measured. 
Even though the observed magnetic moments do not satisfy the SU{3) relations 
very well, a more accurate treatment is unnecessary in view of the other uncer- 
tainties involved. 

The complete determination of the 6j, parameters requires an additional 
assumption to relate the experimentally constrained bi to the unconstrained 6f. 
We impose 6f = bi for i = 1,2 (63, 6| do not enter the EDM expression). This 
choice seems natural starting from the Lagrangian (^) and is our second main 
assumption. 

The linear combinations Bobirris and i?o&2'^s enter the mass terms of the 
baryons as implied by (p. Defining rh = [rriu + mdj/^ we use the combinations 
of baryon masses which are not affected by the isospin breaking effect at tree 
level. They are 

rriN = = m - ABobim - 45o62m, (= 938.91897(28) MeV) 

Mso = ^ ^ — — =m- 45o(6i + b2)Th (1192.55(8) = 1193.41(5) MeV) 
Mh = + Ms- ^ ^ _ 4^^5^^^ _ 4Bob2rh (= 1318.07(11)) MeV^) 

Ma = m--Boibi + b2)im + 2ms) (= 1115.57(6) MeV"), (14) 
3 

where m = itlb — 4i?o&3(2m + m^) takes into account the contribution from the 
63 term in (^). The values in brackets are the latest experimental determinations 
0]. Using the experimental values of the four masses mjv, M^o, M=, Ma in (|T^) 
we can determine the combinations B^bi and i?o&2 with fixed fh and (or 
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alternatively the combinations Bobirris and Bob2ms if we approximate m = 0). 
Using the set (M^o, rriN, Ms) we get 

M^o — miy = ABob2{ms — fh) 253.63 
Mso - Ms = 4So6i(m, - m) ~ -125.52 

(15) 

and with m = 6 MeV, m, = 175 MeV we get 2Bobi = -0.3714 and 2^062 = 
0.7504. These values give M\ ~ 1107.15 MeV, with m extracted from M^o, 
which is a reasonable approximation of the real value. Alternatively if we use the 
set (Mso,mjv,MA) the numbers change to 2Bobi = -0.4088 and 2^062 = 0.7504 
and a slightly too small value for the S mass Ms = 1054.38. 

We still need an additional constraint to fix bi, 62 together with the coefficients 
of the weak Lagrangian 7, 6, or equivalently a and b. The latter enter the weak 
non leptonic hyperon decay amplitudes. There are seven measurable amplitudes: 

— >■ riTT^, S+ prr^, A mr^, A pvr^, —>■ Att^, 2° Att" and three 
isospin relations both for S-wave and P-wave amplitudes. For the chosen four 
independent S-wave amplitudes at tree level one has (for the standard definition 
of the S and P-wave amplitudes in weak hyperon decays see e.g. [0): 

v4(^)(S+ ^ n7r+) = (0.06 ±0.01) 

^nvr") = (1.88 ±0.01) 

A^'\A^pn~) = ^ (1.42 ±0.01) 

A^'\E-^An-) = (-1.98 ±0.01), (16) 

where the last number in parenthesis on the r.h.s. is the corresponding experimen- 
tal value in units of GFm'^+ (this is in agreement with ref. [|1^] since experimental 
values for the decay parameters of hyperon non leptonic decays are unchanged 
since the Particle Data Book of 1990 0). We use the values a = -0.58 ± 0.21 
and b = 1.40 ± 0.12 in units of Gpml+f^ (/^ ~ 132 MeV) [0 determined with 



a tree level least squares fit of the seven measured S-wave amplitudes. These 
determine the combinations 

25o76i = a + b = 0.82 Gprnl+f^ 

2BoSb2 = a-b = -1.98 a pml+U, (17) 

with Gp^'^+fTT = 2.98 ■ 10^^ GeV. Using instead only the measured values of the 
decay amplitudes of the two processes mx~ and A pir^ for a tree level 

determination of a and b one gets: a = —0.54 and b = 1.34. 
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The combinations ( |T7| ) determine the values of the parameters At defined in 
([n]). We obtain 

Aj: = -V2Bo6b2^lAOGFml+U 

Aa = ^Boi5b2-2-fbi)c^-lA8GFml+U, (18) 

while the numerical values obtained by using directly the experimental values 
for — s> niT^ and A — * pir^ decays are 1.33 Gpfn^+f-n and —1.42 Gf^^+Zw 
respectively. Using then hi = h\ and the values in (T3) we also get 



A%, = 2V-2B,l,l^ ^^f-'"" ^1.06 

y lynis — m) 




A\ = -2W^B„(6=-26;) = -?%i^5^.-1.22, (19) 

V o V o(wis — mj 

while for A\ one gets —1.28 if using the experimental values of m^v, and mA. 

The last parameter to be estimated is a— a*. In terms of the Wilson coefficient 
function of the effective four-quark operator we have 

= -HdRdL) ^ Im Ce, (20) 

where 

CM = V^,V:,zM - VtdV*yM (21) 

and 

Im Cg(/i) = -Im VtdV* yM = C23S23S13 sin(5i3 yein)- (22) 

The estimate of the size of the coefficient ye^fi) is affected by large uncertainty. 
We use the renormalization scheme independent definition in [|l2l where the next- 
to-leading corrections at a given /i to the effective hamiltonian are shifted into 
the Wilson coefficient functions. As noticed there the coefficient yei^fj.) is a very 
sensitive function of Aj^jg and being next-to-leading corrections sizable. We use 
the approximate value yQ ~ —0.13 at ~ 1 GeV and A^jg ~ 300 MeV. The CKM 
matrix elements are |T^ S23 = |Kb| = 0.040±0.005, S13 = \Vub\ — 0.0032 extracted 



from iV^fe/Kbl = 0.08±0.02. We approximate cosines to unity and put sin^is ~ 1 
0. Using for the scalar quark condensate {dudi) = — 1/2 ■ (0.235)^ GeV^, this 
gives (a - a*)/2 = -i 5.48 x lO'^^ Q^y 

For the final prediction we use As = 1.40 Gf^^^+A, Aa = —1.48 Gprnl+fT,, 
Gp^'i+fn = 3.0 X 10^^ GeV and the experimental values for the magnetic mo- 
ments and baryon masses in the EDM formulas (^ and ([T^) . We use /^aso = 
-1-1.61. We predict the following value for the neutron EDM: 

C = - "-"7 _9.i8 X 10-s GeV- 
2mN 2 



8 



+7.89 X 10-^ GeV-^ 
-34.02 X 10"^ GeV-^ 
-20.04 X 10"^ GeV-^ 
-20.19 X 10~^ GeV'^ 

-15.87 X 10-^ GeV-^Y (23) 

where the numerical value in each line corresponds to the relative expression in 
(^. This shows the absence of large cancellations. This gives 

^ 5.3 X 10~^2 e ■ cm (24) 

For the proton we have: 

dP ^ -3.6 X 10"^^ e ■ cm, (25) 

where we used the experimental values for A/i2+ = 1.458 in units of nuclear 
magnetons, Mj^+ = 1189.37(6) MeV. Both the neutron and proton electric dipole 
moments acquire the opposite sign if we use instead 6f = — foj. 



4 Power counting and loops 

The purpose of this section is to derive the power counting rules for quantum 
corrections to the tree level electric dipole moment contributions. They can 
be consistently derived within the Heavy Baryon chiral perturbation expansion 
(HBChPt). 

One comment is in order concerning the heavy baryon mass limit of our tree 
level contribution to the electric dipole moment term. In the ordinary Heavy 
Baryon ChPt the EDM term is one of the possible counterterms which appears 
in the tree level Lagrangian at order in the derivative expansion. As an example 
the first tree level diagram shown in Fig.[^ of the full baryon theory leads to the 
following contributions in the heavy baryon mass limit 

-2i^ 2 2 + 2«r2 V- 26 

- ^pole [mi - m^^ie) [mi - mpole) 

These terms are enhanced by one inverse power of the baryon mass splitting 
respect to the ordinary EDM counterterm appearing at order in the derivative 
expansion. 

The heavy baryon mass limit in ( pU]) of the tree level contribution to the 
electric dipole moment derived in the full theory shows that the leading term 
appears at order p in the derivative expansion, while the usual first counterterm 
to the electric dipole moment appears at order p^ both in the full theory and 
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(a) (P) (y) 



Figure 2: One loop diagrams of class a) which contribute to the neutron and 
proton electric dipole moments. The circle vertex is the parity conserving (p.c.) 
vertex. The box is the parity violating (p. v.) vertex. The photon insertion in 
and (7) is the magnetic moment vertex. Internal baryon propagators can be 
also decuplet. 

in the HBChPt. This power counting for (p6|) is consistent with the fact that 
baryon propagators count as 1/p in the derivative expansion and that the leading 
parity violating vertex in the HBChPt appears at order p i.e. B^'^^'^^D^By — 
—2B^S^D^Bv. In this case the octet mass splitting, proportional to the off- 
shellness of the baryon propagator, counts as 0{p), while in the strange quark 
mass expansion it is 0{ms). 

The full baryon Lagrangian contributing to the one loop corrections to the 
electric dipole moment term includes: the Lagrangian which mediates \AS\ = 1 
weak interactions, given by (^) with the inclusion of the meson field through the 
substitution Ae —>■ C^Ae^, the usual strong interaction Lagrangian with the inclu- 
sion of meson interactions starting at order p, the magnetic moment term at order 
p^, the electric dipole moment counterterm at order p"^ and the Lagrangian (|^) 
with the inclusion of the meson field through the substitution A_ ^A_^, A+ — > 
^"''A+^^ Both octet and decuplet states can contribute inside the loop. 

One loop diagrams contributing to the EDM term can be divided into four 
classes: a) corrections to the one loop contribution to the magnetic moment 
through the insertion of the p. v. and p.c. vertices of (^, b) one loop correction 
to the p.c. vertex and to the p. v. vertex in (|^), c) one loop corrections to the 
EDM vertex appearing at order p^, d) one-loop diagrams with the meson loop 
bridging several of the p. v., p.c. and magnetic moment vertices, including the 
case where the p. v. and/or the p.c. vertices emit the meson line. Diagrams with 
the insertion of a photon-meson-baryon-baryon vertex do not contribute. The 
first class is shown in Fig.|^ and the one loop corrections to the magnetic moment 
of baryons in the HBChPt have been derived in |jl3|. 



The magnetic moment one loop contribution is of order p^ in diagram (a) 
and of order p'^ in diagrams {(3) and (7). So the full contribution to the nucleon 
EDM is of order p"^ in diagram (a) and of order p^ in diagrams {(3) and (7). In 
the strange quark mass expansion the counting is somewhat anomalous because 
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Figure 3: One loop diagrams of class c) which contribute to the neutron and 
proton electric dipole moments. The photon insertion is always the EDM coun- 
terterm that appears at order in the chiral perturbation expansion 

the leading tree level contribution starts at order l/m^. One loop diagrams give 
non analytic corrections in the strange quark mass. The one loop in diagram (a) 
gives a correction y^mj to the tree level diagram, while the one loop in diagrams 
{(3) and (7) gives a m^lnm,, correction. 

Diagrams of class b) have the same counting as the corresponding diagrams 
of class a). They include also the electric charge vertex insertion in the proton 
case. Diagrams of class c) start at order and are shown in Fig.§. Diagrams of 
class d) have the same counting as those of a) and only appear at order p^. 

This shows that within the chiral perturbation expansion the tree level contri- 
butions to the electric dipole moment are in fact the leading contributions. One 
loop corrections are suppressed both in the derivative and strange quark mass 
expansion. 

5 Conclusions 

In this letter we have provided a new way of deriving the long distance weak 
contribution to the proton and neutron electric dipole moments. Our final results 
are 

dl ^ ±5.3 X 10^^^ e ■ cm 

dP ^ ^3.6 X 10"^^ e ■ cm . (27) 

The opposite sign is a consequence of the opposite sign in the relevant anomalous 
magnetic moments. These numbers are quite comparable to those derived earlier 
in the quark model and again show that the weak contribution to the electric 
dipole moments is small and of order 10~^^e ■ cm. There are relatively few cancel- 
lations involved in this calculation. We therefore do not expect very large higher 
order corrections. The main uncertainties in the result come from the underly- 
ing assumptions: parity violating terms in the weak \AS\ = 1 Lagrangian are 
negligible, parity violating terms in the strong light quark mass sector are of the 
same size as the parity conserving ones. The other source of uncertainty is the 
estimate of the parameter yelfi) in the Wilson coefficient function of the effective 
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four-quark operator Qq, which we determined according to the renormahzation 
scheme independent definition in [|^. This result should not be added to those 
obtained in the quark model. The value of (qq) is related to the production of 
the constituent quark mass and the contribution as estimated here is thus related 
to the one obtained in the quark model. 
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